Collapse and stable self-trapping for Bose-Einstein condensates with 1/r 6 type 
attractive interatomic interaction potential 
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We consider dynamics of Bose-Einstein condensates with long-range attractive interaction propor- 
tional to l/r b and arbitrary angular dependence. It is shown exactly that collapse of Bose-Einstein 
condensate without contact interactions is possible only for b > 2. Case b = 2 is critical and re- 
quires number of particles to exceed critical value to allow collapse. Critical collapse in that case is 
strong one trapping into collapsing region a finite number of particles. Case b > 2 is supercritical 
with expected weak collapse which traps rapidly decreasing number of particles during approach to 
collapse. For b < 2 singularity at r = is not strong enough to allow collapse but attractive l/r b 
interaction admits stable self-trapping even in absence of external trapping potential. 

PACS numbers: 03.75.Kk, 03.75.Lm 



The dynamics of Bose-Einstein condensate (BEC) with 
short-range s-wave interaction have been the subject of 
extensive research in recent years [l|, [2|. Condensates 
with a positive scattering length have a repulsive (defo- 
cusing) nonlinearity which stabilizes the condensate with 
the help of external trap. Condensates with a negative 
scattering length have an attractive (focusing) nonlin- 
earity which formally admits solitons. However, with- 
out trap these solitons are unstable and their pertur- 
bation leads either to collapse of condensate or conden- 
sate expansion. External trap prevents expansion of con- 
densate and makes solitons metastable for a sufficiently 
small number of atoms. Otherwise, for larger number 
of atoms, the focusing nonlinearity results in collapse of 
solitons. The effect of a long-range dipolar interaction 
on BEC was first studied theoretically 0-0| and more 
recently observed experimentally [8l-ll0j (see also 111. Il2| 
for review). In particular, collapse of BEC with dominant 
dipole-dipole forces predicted based on approximate vari- 
ational estimate [H[ and obtained based on exact analysis 
Q was recently observed in experiment (l3| . 

Here we look for possibility of collapse of BEC due to 
long-range attraction vs. formation of stable self-trapped 
condensate for a general type of long-range interaction 



V(r) 



/(n) 



b > 0, 



(1) 



where /(n) is an arbitrary bounded function |/(n)| < oo 
and r = (xi, x-i, x$). We do not require /(n) to be sign- 
definite. By attractive interaction we mean that /(n) 
is negative at least for some nonzero range of angles so 
that one can choose a wave function to provide negative 
contribution to energy functional. 

Possible experimental realization of (JT|) are numerous. 
E.g., recent experimental advances allow to study inter- 
action of ultracold Rydberg atoms with principle quan- 
tum number about 100 (see e.g. [T3, (HI)- These interac- 
tions between atoms in highly excited Rydberg levels are 
long-range and dominated by dipole-dipole-type forces. 



Strength of interaction between Rb atoms is about 10 12 
times stronger (at typical distance ~ 10/zm) than inter- 
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action between Rb atoms in ground state (see e.g 
for review). Strength and angular dependence of inter- 
action between Rydberg atoms can be tuned in a wide 
range [ljl E.g., spatial dependence for Rb with 

principle quantum number about 100 can be oc 1/r 3 for 
r < 9.5/im and oc 1/r 6 (van der Waals character) for 



> 



9.5/im 15j. Short-range s-wave scattering interac- 



tion is limited to much smaller distance ~ few nm so that 
the range of dominance of long-range interaction poten- 
tial is quite high. Another possible form of long-range 
attractive interaction is gravity-like 1/r potential which 
is proposed to be realized in a system of atoms with laser 
induced dipoles such that an arrangement of several laser 
fields causes cancelation of anisotropic terms 17| . Terms 
oc 1 /r 2 are also possible [17| . 

The mean field BEC dynamics is governed by a nonlo- 
cal Gross-Pitaevskii equation (NGPE) 
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where $ is the condensate wave function, the contact 
interaction is oc g — 4iTh 2 a/m, a is the s— wave scattering 
length, m is the atomic mass, ujq is the external trap 
frequency in the x\ — xi plane, 7 is the anisotropy factor 
of the trap, and the wavefunction is normalized to the 
number of atoms, / d 3 r |'i , | 2 = N. Contact interaction 
term can be also included into potential V(r) as §<5(r) 
but we have not done that because we focus here on effect 
of long-rage potential (JXJ) . If V(r) = then a standard 
Gross-Pitaevski equation (CPE) [l| is recovered. 



NGPE ([2]) can be written through variation iK— - 



j^r of the energy functional 



at 



E = E K + E P + E NL + E R , 



(3) 



2 



which is an integral of motion: ^j- = 0, and 
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K 



2m 



|V*| 2 d 3 r, E NL = ^ / |*| 4 d 3 r, 



Ep = J -^mw 2 (x 2 + x 2 +j 2 x 2 )\^\ 2 d 3 r, 
Er = ^J |*(r)| 2 y(r-r')|*(r')| 2 d 3 rd 3 r'. 



(4) 



Consider time evolution of the mean square radius of 
the wave function, (r 2 ) = / r 2 \^\ 2 d 3 r/N. Using in- 
tegrating by parts, and taking into account vanishing 
boundary conditions at infinity one obtains 



d t (r 2 ) 



2mN 



j 2i.r^i), v|'' - $*d Xj V)d 3 r, (5) 



where dt = Mr, d Xi = and repeated index j means 



summation over all space coordinates, j = 1, . . . , 3. After 
a second differentiation over t, one gets @ 



d 2 (r 2 ) = 



1 



2mN 



8E K - SE P + 12E NL 



-2 J |1r(r| 2 |*(rf (^ck, + ^ )U(r - r')d 3 rj , (6) 
which is called by a virial theorem [6| similar to GPE 

U-Hi. 

It follows from (JTJ) that (xjd Xj + x'fi^V^ - r') = 
-bV(r — r') and using ((3]) we rewrite ([6]) as follows 
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AbE + (8 - 46)^ - (4 + 26)mw 2 7V (r 2 ) 
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-(4 + 26)mw 2 iV( 7 2 - l)(x 2 ) + (12 - Ab)E NL .(7) 



Here both the local nonlinear term Enl and the nonlocal 
nonlinear term Er are included into the energy E. Catas- 
trophic collapse of BEC in terms of NGPE means a singu- 
larity formation, max 1^1 — > oo, in a finite time. Because 
of conservation of N, the typical size of atomic cloud near 
singularity must vanish. The virial theorem ([7]) describes 
collapse when the positive-definite quantity (r 2 ) becomes 
negative in finite time implying max I'J'I — > oo before (r 2 ) 
turns negative. The kinetic energy Ek diverges at col- 
lapse time which follows from divergence of the potential 
energy for max | VP | — > oo together with conservation of 
the energy functional E. Another way to see divergence 
of Ek is from uncertainty relation Ek > y— (9/4) N/ (r 2 ) 
(see [H,[2(|) for (r 2 ) — > 0. Generally (r 2 ) may not vanish 
at collapse (e.g. if there are nonzero values of \^f\ away 
from collapse center) but Ek diverges at collapse time for 
sure because of max 1^1 — s> oo. We use below divergence 
of Ek as necessary and sufficient condition of collapse 
formation while vanishing of (r 2 ) is only sufficient condi- 
tion for collapse. 

NGPE is not applicable near singularity and another 
physical mechanisms are important such as inelastic two- 



and three-body collisions which can cause a loss of atoms 
from the condensate [lj . In addition, multipole expansion 
used for derivation of the dipole-dipole-type potential is 
not applicable on a very short distances (few a few Bohr 
radii). However, as explained above, NGPE with po- 
tential (fTJ) is a good approximation for a wide range of 
typical interatomic distances. 

Consider case 2 < b < 3. Then one immediately ob- 
tains from equation §7$) that d 2 (r 2 ) < Integrat- 
ing that differential inequality over time we get that 



(r 2 ) < 



3E +2 
mN 



t Z + 



dt{r 2 )\t=ot + (r 2 )| t=0 - If E < we con- 
clude that (r 2 ) — > for large enough t which provides a 
sufficient criterion of collapse of BEC. Condition E < 
is sufficient but not necessary for collapse. Using gener- 
alized uncertainty relations between Ek , N,(r 2 ), d t {r 2 ) 
0, 20;] one can obtain much stricter condition of collapse 
which is outside the scope of this Letter. 

Below we assume g = 0. Choosing e.g. initial condi- 
tion as ^| t =o = ^e^/^i) gives E = _*£n + 

7T' 1 / 2 2 4 - b fN 2 r- b T(3°/2 - 6/2) for lj = and /(n) = 
Const = f. If the constant r is small and b > 2 or if 
N > 3h 2 /(l6mf) and 6 = 2 then E < 0. It means that 
for 2 < 6 < 3 long-distance potential alone is enough 
to achieve catastrophic collapse of BEC. Another par- 
ticular example was considered in Ref. @ for the case 
of dipole-dipole interaction potential with all dipoles ori- 
ented in one fixed direction. Note also that b = 3 is a 
border between short-range potentials (for 6 > 3) and 
long-range potentials (for 6 < 3) in 3D (one needs con- 
vergence of J r<r V(r)d 3 r to have short-range potential, 
where r c is a cutoff at small distances). Case 6 = 3 also 
requires that integral of V(r) over angles to be zero to 
ensure convergence of (p} at small r (as is the case |6[ for 
dipole-dipole interactions with fixed direction) otherwise 
we would have to introduce cutoff at small distances and 
potential would loose general form (fTJ). 

Case 6 > 3 generally requires introduction of cutoff 
at small scales. Because J r<r V(r)d 3 r is finite in that 
case we generally have situation very similar to standard 
^-correlated potential [l[ . 

Now we prove that for 6 < 2 collapse is impossible 
for .g = because singularity of (fTJ is not strong enough. 

We use the inequality / j|^d 3 r < 4 J" |W(r)| 2 ci 3 r 0, 
which holds for any r'. Using now Holder's inequality we 
generalize that inequality (assuming 6 < 2) as follows 
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Using now boundness of / : /(n) < f m = max|/(n)| in 
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equation 



FIG. 1: Schematic of the function P(E K ) defined in (TTOJ) 
(solid curve). Equation P(Ek) — E (E corresponds to dashed 
line) has either one of two roots for Ek > depending on sign 
of E. Eg is the largest root. 



(JTJ and inequality (j8]) we obtain a bound for Er in (j4j 



Er > -/ m 2 b - 1 A^ 2 -5 I -E. 



2m \ 2 



J K 



(9) 



which gives a respective bound of E in Q (recall that 
we assume # = 0): 

6 

E > E K - f m 2 b - 1 iV 2 -i ( 2 -^Ek] 2 = P(E K ). (10) 



A function P(Ek) in (I10[) has a minimum for = 

Ef = 2- 2 [/ m 6] 2 /( 2 - b )iV^(2TO/?i 2 ) b /( 2 - b ) resulting in 
a lower bound 



£ > —2- 2 [f m b]—>N—» [-^ 



(11) 



Boundness of the energy functional E from below en- 
sures that collapse is impossible for b < 2. To prove 
that we show that Ek is bounded while collapse requires 
Ek ™> oo. We choose any value of E which satisfy ([TT]) . 
Fig. [T] shows schematically the function P(Ek) from 
(|10[) . Inequality (|10[) requires that Ek < E^(E), where 
Eft{E) is the largest root or equation P{Ek) = E. It 
proves that Ek is bounded for fixed N which completes 
the proof of absence of collapse for b < 2. Particular 
version of that result for 6 = 1 and /(n) = const was 
first obtained in Ref. [25|. Nonexistence of collapse for 
nonsingular potential V(r) was shown previously based 
on approximate analysis in Ref. [261 ] - Proof of nonex- 
istence of collapse for particular example of nonsingular 
potentials with positive-definite bounded Fourier trans- 
form was given in Ref. 27]. These results can be easily 
generalized for any bounded potential similar to above 
analysis. Thus collapse can occur for singular potential 
only and singularity should be strong enough, i.e. b > 2. 

We now look for soliton solution of NGPE ([2]) as 
^(r, t) — A(r)e~ ltlt / h , where a is the chemical poten- 
tial. In that case NGPE ([5J reduces to time-independent 
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A(t) = 0, (12) 



where we again assume g = although generalization 
to g ^ case is straightforward. Equation (fT2"j) is the 
stationary point of the energy functional E for a fixed 
number of particles: S(E — uN) — 0. Multiplying equa- 
tion (fT2|) by A and Xjd Xj A and integrating by parts one 
obtains using (p} and ((SJ that 



Ek,s 

E s 



4-6 



^P,s, E Rs = /J,N S 



4-6' 
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-uN s - — -+2%, 

4 — 6 



(13) 



where subscript "s" means values of all integrals on soli- 
ton solution. Especially simple and interesting is the case 
of self-trapping (ojq — 0) when condensate is in steady 
state without any external trap. All integrals in that 
case depend on number of particles N s only. 

Assume radial symmetry /(n) = Const < in 
(P). Ground state soliton is determined from condi- 
tion that A(r) never crosses zero [H, [29[. To prove 
ground soliton stability we show that it realizes a min- 
imum of the Hamiltonian for a fixed N s . One can 
make inequality ((SJ sharper by minimizing a functional 

J-(^) = A^-i (*2-E K ) V/{^pd 3 r- That minimum 

is achieved at one of stationary points = and af- 
ter simple rescaling one can see that these points cor- 
responds to soliton solutions of the time-independent 
NGPE (fT2| . Among these stationary points the mini- 
mum is achieved at ground state soliton ground- It 
gives a bound J r (^ r ) > minJ r (5') = J- s , ground) which 
is sharper than the inequality ([8]). Most important is 
that following analysis similar to equations (f9"T)- (fTTj) we 
obtain that for any 



E > min E = E. 



s.groundi 



(14) 



i.e. the ground state soliton solution attains the mini- 
mum of E for fixed N . It proves exactly the stability 
of soliton for /(n) = Const. Similar ideas were used in 
a nonlinear Schrodinger equation [28, 29]. Ground state 
soliton was also found numerically for 6 = 1 [3(J . 

For more general /(n) 7^ Const minimum of E is still 
negative if /(n) is negative for a nonzero range of values 
of n. So in that case we expect that ground state soliton 
solution attains that minimum and, respectively is stable. 
If /(n) > for any n then mini? = 0. It corresponds to 
unbounded spatial spreading of NGPE solution for any 
initial conditions. Self-trapping is impossible in that case 
and solitons are possible for cjq ^ only. 
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FIG. 2: Schematic of E(a) from dT5j) for b > 2, b = 2 and 
b < 2 (see text for details). 



Case b = 2 is on the boundary between bounded and 
unbounded energy functional as can be seen from in- 
equalities (JTUJ) and ([13]), If iV > N S)groun d then E 1 is 
unbounded. If N < N^ ground then E > E Siground = as 
follows from ([13"]) for 6 = 2. Thus N s ^ groun d is the critical 
number of particles for collapse. That critical number is 
similar to collapse of GPE in dimension 2 (2D) (as well 
as similar to critical power in nonlinear optics) |18| . It 
is important to distinguish that critical number from the 
critical number of particles of 3D GPE with lu ^ (HQ. 

To qualitatively distinguish different regimes of col- 
lapse and solitons one can consider, in addition to the 
exact analysis above, a scaling transformations [3l| which 
conserves the number of particles 'F(r) — > a~ 3 / 2l $(r/a). 
Under this transformation the energy functional E (for 
loq = 0) depends on the parameter a as follows 

E(a) = ci~ 2 Ek + a~ b Eji. (15) 

The virial theorem ([7]) and the relations (TT3l) , (ITS"]) have 
striking similarities with GPE if we replace b by the spa- 
tial dimension D in GPE. That analogy suggests to refer 
the case b = 2 as the critical NGPE and b > 2 as the 
supercritical NGPE. Fig. [2] shows typical dependence of 
(TTSl) on a for b > 2, b = 2 and b < 2 assuming Er < 0. 
For b > 2 there is a maximum of E (curve 1 in Fig. [2]) 
corresponding to unstable soliton. Solution of NGPE ei- 
ther collapse or expand. For b = 2 there is no extremum 
and collapse is impossible for N < N s _ groun d (curve 2 in 
Fig. [2]) while condensate can collapse for N > N S:groun d 
(curve 3 in Fig. [2]). Ground state soliton corresponds to 
N = N s groun d and E = locating exactly at the bound- 
ary between collapsing and noncollapsing regimes. For 
b < 2 there is a minimum which corresponds to stable 
ground state soliton (curve 4 in Fig. [5]). 

Solutions of both GPE and NGPE with w = near 
collapse typically consist of background of nearly lin- 
ear waves and a central collapsing self-similar nonlinear 
core. The scaling (fT5|) describes the dynamics of the 
core with time-dependent a such that a — > near col- 
lapse. Waves have negligible potential energy but carry 



a positive kinetic energy E waves ~ E K , waves . The to- 
tal energy E = E coUapse + E waves is constant, where 
Ecoiiapse is the core energy. It follows from (|15j) that 
for b — 2 one can simultaneously allow conservation of 
N and E co u apse so that negligible number of waves are 
emitted from the core. This scenario is called a strong 
collapse. If b > 2 then the term cx a~ b in (p~5|) dominates 
with E co u apse — s> — oo as a — > 0. Then the only way to 
ensure conservation of E is to assume strong emission of 
linear waves. Near collapse time only vanishing number 
of particles remains in a core (of course all that is true 
until NGPE losses its applicability) which is called weak 
collapse [3l|. Strong collapse was shown to be unstable 
for supercritical GPE [3l|. It suggests that collapse for 
b > 2 is of weak type. 
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